Exact coherent matter-wave solitons induced and controlled by laser field 
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We find a set of exact solutions of coherent bright solitons in the quasi-one-dimensional (ID) 
Bose-Einstein condensate (BEC) trapped in a harmonic potential, by using a Gaussian laser well 
(barrier) with oscillating position to balance the repulsive (attractive) interatomic interaction. The 
bright solitons do not deform in propagation and are controlled accurately by the laser driving which 
resonates with the trapping potential. The solitonic motion is more stable for the repulsive BEC 
than that of the attractive BEC. The results reveal a different kind of soliton trains compared to 
that reported recently in Phys. Rev. Lett. 100, 164102 (2008) and suggest an experimental scheme 
for generating and controlling the coherent matter-wave solitons. 
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Soliton in a Bose-Einstein condensate (BEC) is a kind 
of important nonlinear phenomena, which has been in- 
vestigated experimentally and theoretically The 
formation and propagation of BEC solitons were ob- 
served by magnetically tuning the atom- atom interaction 
from repulsive to attractive P, Q • The theoretical works 
demonstrated that the matter-wave bright solitons can 
be created in a BEC, through the modulational insta- 
bility 0, 0| and quantum phase fluctuations Prop- 
agation feature of the solitons is the breathing oscilla- 
tion which is mostly controlled by the harmonic trap 
fisj - In the harmonic trapping case, the Gaussian 
ansatz was used to fit the profiles of solitons [1, 0, [3, [3 ■ 
The Gaussian-shaped optical potentials have been ap- 
plied for investigating the BEC solitons [l^, [13, [3 ^^'^ 
quantum tunneling [l^, (23, S3- worth noting 

that the balance between nonlinearity and dispersion was 
found in the seminal study of soliton . Recently, the 
new balances between the atom-atom interaction and the 
Gaussian and/or periodical potentials are demonstrated 
[TtI [2^ 1 . For some special forms of external potential and 
interaction intensity, the exact soliton solutions in BECs 
have also been reported [H, [H, [H, 113 . 

The quantum states governed by the linear Schrodinger 
equation with inseparable space-time variables are very 
important but had to find. The coherent state of a har- 
monic oscillator is a nice example of such states, which 
has been widely applied to physics and optics [H, [2§| 
and is also extended to the case of wavepacket trains 
[30| . Can the coherent wavepacket trains exist in 
a harmonically trapped BEC system governed by the 
Gross-Pitaevskii equation (CPE)? This is usually impos- 
sible, because of the nonlinearity in CPE. However, when 
we employ the laser field to balance the nonlinear term, 
seeking exact coherent states of the GPE could become 
possible. Demonstrating the exact coherent state of GPE 
and its experimental feasibility is our main motivation in 
this paper. 



By using the balance technique and applying the os- 
cillating Gaussian lasers, we find n exact solutions of co- 
herent soliton trains in the quasi-lD BEC. When n = 
is considered, for the attractive or repulsive BEC the in- 
teratomic interaction is balanced by the Gaussian barrier 
or well respectively, and in the both cases the soliton so- 
lutions possesses the same form of coherent state. The 
coherent bright soliton oscillates like a classical harmonic 
oscillator with the trapping frequency, which agrees with 
Strecker's experimental results However, compared 
to the deformed soliton trains observed in experiment, 
our bright solitons have different properties, namely their 
shapes are kept in propagation, their behaviors are con- 
trolled accurately by the laser field and their motions 
possess more stability for the repulsive BEC rather than 
the attractive one. Based on the capacity of current ex- 
periments, such bright solitons can be observed in a BEC. 

We consider a BEC consisting of N identical Bose 
atoms and being transferred into a cigar-shaped har- 
monic trap. The potential that takes into account the 
combination of the magnetic trap with the laser sheet 
reads V{x,t) = ^mculx'^ + VL{x,t) with V>,(x, t) contain- 
ing the Gaussian-shaped laser potential [H, [21I, [l^ . Let 
the transverse frequency ujr be much greater than the ax- 
ial frequencies uj^ , the dynamics of the system is governed 
by the quasi-lD GPE 0,[i| 



2m dx^ 



[V{x,t)+g[M']i^, (1) 
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where we have assumed the transverse wave function be- 
ing in ground state of a harmonic oscillator such that 
the quasi-lD interaction intensity related to the s-wave 
scattering length , atomic mass m and number of con- 
densed atoms TV reads [sH — Nmuorgo/ {^irh) — 
iNhi/JrCLs for the normalized wave-function ip. The norm 
IV'P is the probability density and A^|'0P the density 
of atomic number. Setting 1,. — ^Jh/ {mcur), Ix = 
^Jh/{mijJx), we normalize the time, space, wave func- 
tion and laser potential by Ix, \/Qlx and hwx re- 
spectively, then the interaction intensity becomes gid — 



2 



2Nu>ras/ {lUxIx) and the dimensionless GPE reads 



'2 9^ 



-x'' + VL{x,t)+giM'' i^- (2) 



In order to seek the extended coherent state of Eq. (2), 
we have to use the balance condition [23| 

VL{x,t)+giM^ = (3) 
to transfer Eq. (2) to the Unear Schrodinger equation 
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(4) 



where fj, denotes an undetermined constant determined 
by the laser profile and is in units of fuvx- Obviously, 
by the balance condition we mean that the external op- 
tical potential and the internal interaction reach into an 
indifferent equilibrium experimentally, namely their sum 
equates a constant [l^l- The exact solutions of Eq. (2) 
must obey the balance condition (3) and linear equation 
(4) simultaneously. Therefore, only the properties com- 
mon to the nonlinear Eq. (3) and linear Eq. (4) can be 
kept in the balance solution. For the Gaussian potential 
VL{x,t) the soliton solution of the nonlinear Eq. (3) and 
coherent-state solution of the linear Eq. (4) can be in the 
same form and could coexist thereby. 

The exact solutions of extended coherent states of Eq. 
(4) read 

iPn = i?„(x,t)e''^"(^'*\ n = 0,1,2,- •• ; 
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+ fi + n^t + XQXSmt ^ —XQSm2t , (5) 



which can be proved by inserting it directly into Eq. (4). 
Here denotes the Hermite polynomial with xq being 
the amplitude of the center position of Gaussian packet, 
which can be adjusted by the oscillating amplitude of 
laser position. Applying Eq. (5) to Eq. (3) shows the 
profile of required laser field 



VL = ^i- 



7^2 



Clearly, for constant interaction gid and n = {H^ = 1) 
case Eq. (6) describes the oscillating Gaussian potential 
which agrees with that used in a;o = it be- 

comes the well known time- independent form [13, E, l2lj | 
which leads Eq. (5) to the exact stationary state of Eq. 
(2). For any n Eq. (6) denotes a multi-well (-barrier) 
which includes the well-known double- well with n — \. 
The absolute value \gid\ determines the required laser in- 
tensity. Obviously, Eq. (5) describes the oscillating single 
soliton for n = and the soliton trains for n > [l|, [ij] ■ 
The oscillating amplitude Xq and frequency lo = \(ijJx) 
of the solitonic center are controlled by the oscillating 
laser strictly. The required laser frequency ojl is equal 



to the trapping frequency oJx, that means resonance be- 
tween the trapping and driving potentials. The same 
frequency fixes width and height of the Gaus sian laser 
and solitons, through the axial length = \/ h/{muJx) 
of harmonic oscillator. On the other hand, rewriting 
Eq. (5) as ipn — t)e~*^"*/'', we obtain the Flo- 

quet quasienergy En — (| + + adjusted by the laser 
parameter /i. The corresponding Floquet state obeys 

(l)n{x,t + 2tt) = 0„(x,t). 

It is interesting noting that the soliton solutions of Eq. 
(5) possess the same form for the attractive and repulsive 
interatomic interactions. However, for different interac- 
tions the solitons are associated with different shapes of 
the required Gaussian lasers. Taking n = as an exam- 
ple, Eq. (6) exhibits that the laser barriers correspond 
to the attractive interaction with gid < and the laser 
wells are associated with the repulsive one with 171^ > 0. 
When we consider N = 10^ ^Li atoms with mass ni 
being 7 times of the proton mass rrip and take the ex- 
perimental parameters [H = ±1.5nm, lu^ — 20Hz, 
ujr = 800Hz, the harmonic oscillator lengths and the in- 
teraction intensity become Ix = \/40Zr = 21.22/^m and 
gid = 2Nujras/ (lUxIx) — ±56.55. In Figs, la and lb we 
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FIG. 1: Spatiotemporal evolutions of the atomic density plot- 
ted from Eq. (5) for the laser parameter xq ~ 10{lx) = 
212.2/im and (a) n = 0, (b) n = 1. At any time Fig. la 
displays a single bright soliton and Fig. lb exhibits a soliton 
pair, both oscillating their centers. 

show the spatiotemporal evolutions of the atomic den- 
sities i?g(a::,t) and Rf{x,t) respectively for the laser pa- 
rameter a;o — 10(^2;) = 212.2/j,m. The former is a single 
Gaussian wave and the latter is a double Gaussian at any 
time for both gid > and 171^ < cases. These bright 
solitons oscillate their centers with amplitude xq. Here- 
after, all the parameters in any figure are dimensionless. 
The laser potential VL{x,t) of Eq. (6) with parameters 
/i = 10, Xq = 10 is plotted as in Fig. 2 for four sets of 
parameters respectively. In Figs. 2a and 2c with positive 
gid, we observe the single and double wells respectively 
for any time. The wells oscillate their center positions as 
the increase of time. From Figs. 2b and 2d with negative 
did, we find that at any time the laser potentials describe 
the single and double barriers with oscillating centers. 

To see the details of the soliton motions and to ana- 
lyze the stability of the system, in Fig. 3 we show spa- 
tial profiles of the soliton and total potential functions 
V{x, t) = ^x'^ + VLix, t) at several different times for the 
same parameters with Fig. 2a. It is revealed that in the 
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FIG. 2: Spatiotemporal evolutions of the potential func- 
tions plotted from Eq. (6) for (a) n — 0, gid = 56.55, 
(b) n = 0, gid = -56.55, (c) n = 1, gid = 56.55, (d) 
n — 1, gid = —56.55. The potential wells and barriers are 
found respectively for the different cases. 
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FIG. 3: Spatial profiles of the soliton (solid curves) and po- 
tential functions (dashed curves) for the parameters of Fig. 
2a and at (a) t = 0, (h) t = tv/2, [c) t = tv and (d) t = 57r/4. 
It is observed that the soliton keeps its shape and oscillates 
its position. The soliton center is located on the center of a 
potential well for any time. 



time evolution the potential deforms between the single 
well and double well and the density soliton keeps its ini- 
tial shape. When the initial time is taken as to = 0, from 
Eq. (5) we plot the initial profile of density times interac- 
tion intensity gidR^ as in Fig. 3a, which is a bright soli- 
ton centered at x = 10. With time increasing to t = tt/2 
and t — -K, the soliton propagates toward left to posi- 
tions a; = and x = —10 as in Figs. 3b and 3c. In such 
a time interval, the potential deforms from a double well 
to a single well then to another double well. In the next 
deformation period tt < LDxt < lix of the potential, the 
soliton propagates toward right as shown in Fig. 3d and 
will back to the initial place. The propagation property 
is in good agreement with that of Strecker's experiment 
and the potential in Fig. 3a has the profile illustrated in 
previous work [3^. The soliton center is always located 
on the center of a potential well. This implies that at any 



time as a quasi-particle the soliton falls into a potential 
well. According to the well-known criterion of dynamical 
stability, such soliton motion is stable. In the transporta- 
tion process of BEC soliton, the laser potential Vi,(x,t) 
plays a role of optical tweezer [33i, i34] • 
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FIG. 4: Spatial profiles of the soliton (solid curves) and poten- 
tial functions (dashed curves) for the same parameters with 
Fig. 2b and at (a) t = 0, (b) t = 7r/2, (c) i = tt and (d) 
t = 57r/4. At t = and t = tt the soliton is not located on 
the position of minimal potential. For t — 71/2 and t — 5tt/4: 
the soliton is located on the position of potential barrier. 

When the interaction intensity is changed from gid — 
56.55 to g\d — —56.55, Fig. 3 is correspondingly changed 
to Fig. 4. In the latter figure, the solitonic shape and 
evolution have no change. However, the potential de- 
forms with different pattern compared to Fig. 3 such 
that the soliton is no longer located on the center of a 
potential well. Particularly, Figs. 4b and 4d exhibit that 
sometimes as a quasi-particle the soliton lies at the top 
of potential barrier. The corresponding soliton motion 
may be dynamically unstable thereby. 

Similarly, in the double-soliton case with n = 1, by 
comparing Fig. lb with Figs. 2c and 2d, we find that 
the solitons fall on the optical wells or barriers for the 
repulsive or attractive BEC. Therefore, the exact soli- 
ton pair is dynamically stable for the repulsive BEC or 
unstable for the attractive one. 

In conclusion, we have investigated the repulsive and 
attractive quasi-lD BECs held in the combination po- 
tential of the magnetic trap and the Gaussian laser sheet 
with oscillating position. It is demonstrated that when 
the laser potential balances the interatomic interaction, 
the exact bright soliton trains can be generated. The cor- 
responding n soliton solutions agree with the extended 
coherent states of harmonic oscillator. The soliton trains 
fit the periodical motions of laser centers and keep their 
shapes of Gaussian waves, that agree with Strecker's ex- 
periment partly [l|. The solitonic width, height, oscillat- 
ing amplitude and frequency are controlled by the laser 
field accurately. The required optical potentials contain 
the Gaussian wells and barriers for the repulsive and at- 
tractive BECs respectively, which resonate with the trap- 
ping potential. For n = case the optical well is similar 
to the quantum dot generated by a focused beam of red- 
detuned laser light "3311. Although the solitonic profile of 
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fixed n is same for the both interaction cases, the soliton 
of repulsive BEC is more stable than that of the attrac- 
tive BEC. It is worth noting that in some time intervals 
the BEC solitons are transported toward fixed direction 
and the laser field behaves like a quantum tweezer re- 
alized in previous works [s^ [13] ■ Therefore, the exact 
bright solitons could be observed and controlled experi- 
mentally by oscillating the laser position and adjusting 
the system parameters. 

It is also noted that the exact soliton trains are similar 
to the results of Ref. [l^. However, the required external 
potentials and interatomic interactions are different for 



the both cases. The spatiotemporal-dependent interac- 
tion intensity used in [l^ is not required in our systems. 
This could bring convenience to the experimental obser- 
vation of the soliton trains. 
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